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ABSTRACT 



An optical improvement for angular position sensors, which 
may be used to determine the spatial coordinates of a small 
source of light (or other energy) in a 3-dimensional volume. 
Such sensors normally include a linear photosensitive image 
detector such as a photodiode array or a charge-coupled 
device (CCD). An irregular pattern of parallel slits is 
described which increases the amount of hght gathered 
while avoiding the undesirable characteristics of lens optics 
for this apphcation. One optimal type of irregular pattern is 
the uniformly rcdimdant array. A mathematical correlation 
function together with a polynomial interpolation function 
can determine the displacement of the image on the detector 
and thereby the location of the source relative to one angular 
dimension. Given the locations and orientations of several 
sensors in a 3-dimensional coordinate system and given the 
angles measured by each, the location of the point source 
can be computed. 

20 Claims, 11 Drawing Sheets 
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SYSTEM FOR DETERMINATION OF A Proc. of the 31th International Instrumentation 

LOCATION IN THREE DIMENSIONAL Symposium, San Diego, May 5-9, 1991, p. 163-182. 

SPACE Waldean A. Schulz; "METHOD AND APPARATUS FOR 

THREE-DIMENSIONAL NON-CONTACT SHAPE 

This application is a continuation in part of provisional 5 SENSING", U.S. Pat^Na 5 198377^^ ^M^^^^ 

««rr««t;^„ c*, KT^ /in/n^iiTi <(ua a iqot Farhad Daghighian; "OPTICAL POSITION SENSING 

patent application Ser. No. 60/058,371 filed Sep. 9. 1997. DUOLATERAL PHOTOEFFECT DIODES", 

FIELD OF THE INVENTION Sensors, 1994 November p. 31-39. 

Robert P Burton and Ivan E. Sutherland; "TWINKLE BOX 

This invention relates to an improvement in the optics for A THREE-DIMENSIONAL COMPUTER INPUT 

an electro-optical sensor, such as those that are used as part DEVICE". AFIPS Conference Proceedings 43, 1974, 

of a system that determines the location of one or more point Chicago, 111. 

sources of radiation, such as light, in a three-dimensional The complete disclosures of the above H. Fuchs, Sharon S. 

volume. Welch, and Robert P. Burton references arc incorporated 

herein by reference. 

BACKGROUND OF THE INVENTION Such systems typically include multiple angular location 

sensors. Furthermore, each angular location sensor typically 

Vanous methods and systems have been proposed to includes a linear position detector, such as photo-electric 

measure the location of a radiation source in three- position sensitive detector (PSD), a linear photodiode array, 

dimensional space. By using one or two-dimensional mea- or a linear charge-coupled device (CCD). The linear position 

surements of the location of a source of radiation, typically 20 detector may be a photosensitive semiconductor strip (in the 

measurements that have been made from multiple positions, case of a PSD) or a row of many discrete photosensitive 

the location of the light source in the space can be deter- elements called pixels (in the case of a CCD or photodiode 

mined by calculation. One such method employs multiple array). Area (or two-dimensional) position detectors may 

angle-measuring optical sensors, each of which measures also be used, but the present invention is particularly useful 

the location of the light source with respect to one of several 25 for linear (or one-dimensional) photodiode arrays or CCDs, 

angular dimensions. With the multiple sensors situated which detect an image. The determination of a point in space 

within a three dimensional coordinate system at known using a two dimensional deteaor can be considered similar 

locations and orientations with respect to the coordinate to determining that position by means of two linear position 

system, the 3-dimensional (3-D) coordinates of the light detectors that are disposed at right angles to each other, 

source can be computed relative to that coordinate system. 30 Therefore, the following description will concentrate on a 

In these determinations, the light sources are assumed to be single linear imaging detector that uses discrete pixels. It 

tiny with respect to the size of the volume in which they are should be noted that, while a linear detector measures only 

located. They therefore have been considered to be point one angle (for example elevation) to the point source of 

sources that radiate light outward over a relatively wide radiation, an area detector measures two angles (for 

solid angle. In actuality, there is no such thing as a literal 35 example, elevation and azimuth) simultaneously. However, 

point source of light or other radiation. However, for men- two linear image detectors require only 2N pixels, while an 

suration purposes described herein, a very smairiight source area image detector has N-N' pixels, where N is the number 

in a relatively large volume of space is an adequate approxi- of pixel rows and N' is the number of pixel columns. Since 

mation. For brevity, these approximate point sources of n and N' usually exceed 500, the area detector requires the 

radiation will hereinafter be referred to simply as point 40 readout, storage, and processing of at least several hundred 

sources of radiation or light. times as many pixels. 

Many 3-D methods and systems have been described in Each angular sensor also includes an optical system to 

previous Hterature and some have even been used in practice focus the radiation from the point source into a narrow lineal 

to determine the location of a point source of light in a three real image, and to cause that image to cross the linear image 

dimensional volume. Examples of such prior art are found in 45 detector at an approximately right angle. In reference to FIG. 

the following references: 1 (prior art), a cylindrical lens 27 is typically employed to 

H, Fuchs, J. Duran, B. Johnson, and Zvi. M. Kedem; focus light rays 15 from a point source 10 into a linear real 

"ACQUISITION AND MODELING OF HUMAN image 31 that crosses the photosensitive strip 16 of the linear 

BODY FORM DATA", Proc. SPIE, v. 166, 1978, p image detector 14 at approximately a right angle. 

94-102. 50 Standard spherical lens systems in general cannot be used 

Jean-Qaude Rcymond, Jean-Luc Hidalgo; "SYSTEM FOR with a one-dimensional angular sensor for this purpose, 

MONITORING THE MOVEMENTS OF ONE OR because the focused real image of a point source 10 pro- 

MORE POINT SOURCES OF LUMINOUS duccd by a spherical lens is a tiny spot of light. Because this 

RADIATION", U.S. Pat. No. 4.209,254, Jun. 24, 1980. image is a spot, rather than a line, in most cases it will focus 

Y. Yamashita, N. Suzuki, M. Oshima; "THREE- 55 off the photosensitive strip, that is, far to one side of it. This 

DIMENSIONAL STEREOMETRIC MEASUREMENT is because, although the linear detector's row of photosen- 

SYSTEM USING OPTICAL SCANNERS, CYUNDRI- sitive pixels is long, it is only a few microns wide. Therefore, 

CAL LENSES, AND UNE SENSORS", Proc. SPIE, v. a cylindrical lens 27, rather than a spherical lens, has 

361, 1983, p. 67-73. generally been used in such optical sensors to focus the point 

F. Mesqui, F Kaeser, and P. Fischer; "REAL-TIME, NON- 60 light source 10 into a line image 31. Some tiny portion of 

INVASIVE RECORDING AND 3-D DISPLAY OF THE this line image intersects the linear row of pixels at approxi- 

FUNCnONAL MOVEMENTS OF AN ARBITRARY mately a right angle no matter where the point source is 

MANDIBLE POINT', SPIE Biostereometrics 602, 1985, located in a relative large volume. A compound cylindrical 

p 77-84. lens, consisting of several individual positive or negative 

Sharon S. Welch, Kevin J. Shelton, and James 1. Clemmons; 65 focal length lenses, has sometimes been used to mitigate 

"OPTICAL POSITION MEASUREMENT FOR A optical problems. However, for simplicity only a single 

LARGE GAP MAGNETIC SUSPENSION SYSTEM", simple cylindrical lens 27 is iUustratcd in FIG. 1. 
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As shown in FIG. 2, a simple narrow aperture or slit 22 mination of the location of the point light source. Second, for 

may be employed instead of employing a cylindrical lens to a fixed-focus arrangement of a cylindrical lens, the point 

produce the linear real image, in the same way that a pinhole where the line image 31 meets the photosensitive strip does 

can replace a spherical lens in a conventional camera. A maintain a sharp focus as the point source 10 moves 

single, straight, precision slit 22 focuses light from a point 5 away from the optical centcrline 11. This is especially true 

source 10 onto a line image 32 that approaches a perfect line distance 12 increases, because a standard cylindrical 

image. Such a slit might be a long, very narrow rectangular presents a less circular and more elliptical curvature to 

aperture within an opaque mask, for example. Furthermore, ^^^se light rays 15 entering the lens from an angle that is far 

unlike the cylindrical lens shown in FIG. 1, the slit 22 has ^^'""^.^^ °PV^^^ centerlme of the camera. The 

an infinite depth of field. That is, except for diffraction lo ;csuh is an aberatcd rcahmagc 31 which is m p^^^ 

effects,theJimage32isinsha^f^^^^^^^^^ ^^^S^^^ 

distance of the ponit source 10 relative to the bnear dct^^^^^^ CCD with pixels that individually have diminished 

Furthermore, the image wiU be as straight as the sht is, and intensity). A third disadvantage of the lens 21 is that the 

present day technology can inexpensively manufacture sharpness of focus wiU be affected by the distance of point 

masks with very precise optical apertures. So, as the light 15 source 10 from the lens 21, especiaUy for wider lenses that 

source moves along a path 12 parallel to the longitudinal collect more light. This is the same depth-of-field effect as 

axis of the slit, the point of intersection of the line image 32 for standard spherical lenses with larger apertures, 

and the photosensitive strip 16 remains constant. These disadvantages limit the accuracy of measurement. 

Furthermore, the image remains sharp regardless of the Multiple lenses and light stops can improve the image 

distance (range) of hght source 10 from the slit 22. Further, 20 characteristics, but .such improvements are more costly, 

the angular field of view may be easily changed by varying involve more diflBcult manufacmring, and still do not solve 

the distance between the slit 22 and the photosensitive strip the problems well enough for high accuracy measurement. 

16. As the distance decreases, the field of view increases. The nonlinear distortion can be modeled and corrected in the 

These are significant advantages over lens based optics. computation that calculates the three-dimensional location 

Unfortunately, one sigruficant drawback to the slit 22 is its 25 of the light source from the angle measurements provided by 

Umited Ught gathering abihty. This Umitation produces a multiple sensors. However, this correction is complex 

dim line image 32 in comparison to the image 31 produced P''^'^^^^ '^'^ ^° approximations. This is because the 

by a typical cylindrical lens 27. This limitation restricts the ^^^^ distortion correction for any one sensor cannot be made 

„ „^ „,u:^u *«♦ M^ut using only the data from itself alone. The correction depends 

range at which any practical sensor can detect the pomt ueht ^ r , • j i . • ^ r 

° . ^ . . . r * -J . on knowmg both the azimuth and elevation angles of 

source 10, because, as IS well known, the amount 01 mcident 30 . f^. i- i_* . * .t. i- j- ii u . 

, J ' , ' r .1 . . mcidence of the ught with respect lo the cylmdncal lens, but 

hght decreases inversely with square of the mcreasing ^^^^ ^^^^ measures only one angle with respect to its own 

distance. Even at shorter ranges, where the bnghtness of the j^^^ y^^^^ ^ ^ ^ replacement for the cylindrical 

unage is suffiaent for measurement, the light focused from j^^s does avoid the above problems, its limited Ught gath- 

a single slit still presents a poor signal-to-noise ratio (SNR), ering ability, limited range, and its susceptibility, to dust, 

which limits the reliability of the measurement. 35 smudges, and imperfections are major drawbacks. 

Alternatively, a knife-edge (in effect, a "one-sided slit") or Therefore, there is a need for an optical system that captures 

a very wide slit may be substituted for the slit 22. These the advantages of both slits and cylindrical lenses while it 

options generally are worse than a narrow slit, because they avoids their disadvantages. 

flood the photosensitive detector with more ambient light, Multiple apertures have been used in the prior art, but for 

which competes with the real image of the point source 40 different purposes than those of the instant invention. In the 

while providing no additional substantive position informa- prior art, coded apertures (a patterned form of multiple 

tion. apertures) have been used exclusively with area arrays, such 

A second significant drawback to the slit 22 is its suscep- as in connection with the area CCDs that are used in video 
tibility to dust, imperfections, and smudges. For example, an cameras. The application of such coded apertures has pre- 
opaque particle within the slit but near one edge may cause 45 viously been oriented toward the capture and reconstruction 
a shift in the centroid of the image when the shadow of the of complex two-dimensional images for visualization pur- 
particle falls on the linear photosensitive detector. While poses. They have not been used to improve the precision of 
there are ways to detect when this problem arises, it is one dimensional sub-pixel location measurements, 
undesirable and affects accuracy. Historically, coded apertures have been used to collect more 

In reference to both FIGS. 1 and 2, note that a typical 50 of the available radiation in X-ray imaging, because weaker 
cylindrical lens 27 is more efficient than a slit 22 in con- sources of radiation are desirable from the patient's health 
centrating the light that passes through it and is focused onto perspective, and because conventional lenses cannot be used 
a small area on the photosensitive row of pixels 16. This is to focus X-rays. Coded apertures have also been used in 
because the lens is considerably wider, and more light enters connection with two-dimensional astronomical imaging, in 
the lens 27 than the slit 22. Yet this light is concentrated onto 55 which the sources of light are extremely dim. In the prior art, 
roughly the same area as the image formed by the slit. That the apertures that were used for these purposes were usually 
is one major advantage of using a lens. Another advantage plural pinholes arranged in a pseudo-random, two- 
is that the lens 27 is substantially immune to the effects of dimensional array. 

a moderate amount of dust on it, unUke a single slit aperture The following references are cited as examples of the 
22. .60 stale of the prior art being referred to above: 

However, the cylindrical lens 27 has at least three major Kevin Byard; "SYNTHESIS OF BINARY ARRAYS WITH 

drawbacks. First, the linear real image 31 that is formed by PERFECT CORRELATION PROPERTIES CODED 

a cylindrical lens is not an ideal straight line at all As the APERTURE IMAGING", Nuclear Instruments and 

light source moves abng a path 12 parallel to the longim- Methods in Physics Research A336, 1993, p. 262-268. 

dinal axis of the lens, the point of intersection of the line 65 Kevin Byard; "AN OPTIMISED CODED APERTURE 

image 31 on the photosensitive strip 16 moves a small IMAGING SYSTEM", Nuclear Instruments and Method 

distance. This adversely affects the acciu-acy of the deter- in Physics Research A313, 1992, p. 283-289. 
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Walter C Chiou and Richard C Augeri; "SQUARE ANTI- A still further objective of this invention is to provide a 

SYMMETRIC UNIFORMLY REDUITOANT ARRAY design for the optical apparatus of this invention that mini- 

CODED APERTURE IMAGING SYSTEM", U.S. Pat. mizes manufacturing complexity and yet improves the 

No. 5,606,165, Feb. 25, 1997. performance/cost ratio. 

E. E. Fenimore; "CODED APERTURE IMAGING: PRE- 5 StiU another objective of this invention is to increase the 

DICTED PERFORMANCE OF UNIFORMLY REDUN- amount of location information contained in the image, such 

DANT ARRAYS'*, Applied Optics 17, 22, 1978, p. as the provision of redundant data. 

3562-3570. ^ further objective of this invention is to decrease the 

E. E. Fenimore and T M. Cannon; "CODED APERTURE ^^^cts of non-uniformity of re^onsiveness among the pho- 

IMAGING WITH UNIFORMLY REDUNDANT lO tosensitive elements (pixels) of the detector. 

Ro'^^ii^^^^^^ Other^daddition^objectsofthisinve.^^^^ 

FOR DECOD NG A CODE IMAGE PRODUCED BY ' consideration of this specificaUon mcluding 

PENETRAHNG RADIAHON, SUCH AS X-RAYS", ^^^^.^.^^^ ^^^^^^ '''' ^ appended hereto. 

U.S. Pat. No. 4,191.890, Mar. 4, 1980. 15 Satisfying some or all these objectives will d^ectly 

Stephen R. Gottcsman and Edward J. Schneid; "PSEUDO- ^^^'^^^ measurement accuracy and reliabihty of indi- 

NOISE PRODUCT CODED APERTURE ARRAYS vidual optical sensors, and therefore the accuracy of a whole 

AND METHOD FOR DESIGNING SAME", U.S. Pat. three-dimensional system employing them will be 

No. 5,036,546, Jul. 30, 1991. unproved. 

S. R. Gottesman and E. J. Schneid; "PNP— A NEW CLASS 20 invention described herein provides an improvement 

OF CODED APERTURE ARRAYS", rra/isacftofts ^ '^c optics used in conjunction with an electro-optical 

on Nuclear Science 33.1, 1986 February, p. 745-749. ^^^r image detector. Together, these constitute a sensor that 

Philippe Laudel and Jean-Pierre Roques; "RESOLUTION measures the one-dimensional angular location of a light 

OF STRONG SOURCES FOR A GAMMA-RAY TELE- source. This improvement is also intended to be applied to 

SCOPE USING CODED APERTURE IMAGING", 25 *h® individual sensors of more complex, multi-sensor sys- 

AppUed Optics 27:20, Oct, 15, 1988 p. 4226-^230. ^^ms. These more complex systems may employ several 

Masaru Matsuoka and Yoshili Kohmura; "A NEW CON- such sensors to measure the angular location of a point light 

CEPT OF X-RAY MICROSCOPES WITH A CODED source from several locations and from these data to calcu- 

APERTURE IMAGING MASK", Jpn, J. AppL Phys, late the coordinates of the Ught source within a three- 

MJpn. J. AppL Phys, 34:1:1, p. 372-373. 30 dimensional coordinate system. Simpler versions of the 

Jean in 't Zand; "CODED APERTURE CAMERA IMAG- apparatus of this instant invention are adapted to make 

ING CONCEPT", Internet WWW site http:// measurements within a one- or two-dimensional coordinate 

lhcawww.gsfc.nasa.gov/docs/cai/coded_intr.html, Oct. 7, system. Furthennore, versions of the apparatus of this inven- 

3^997 tion might employ multiple light sources attached to a rigid 

Included in some of the above references is the description 35 body so the location and orientation of the body as a whole 

of a particular form of coded aperture called a uniformly ^ig^^ ^e determined. 

redundant array (URA). While the references employ two- To achieve these and other features, the invention pro- 
dimensional URAs for imaging purposes, the present inven- vides an improved optical apparatus for receiving radiation, 
tion uses one-dimensional URAs specifically for angular preferably light or electromagnetic radiation, from a point 
position sensors. URAs have particular properties, which 40 radiation source and determining the angular location of the 
improve the results of applying a mathematical correlation source about an axis of the apparatus. The apparatus com- 
to the image in order to determine its exact location. Unlike prises a linear imaging detector that itself comprises an array 
many arbitrary aperture patterns, which may generate mul- of many light sensitive elements (pixels) that are adapted to 
tiple local peaks in the correlation, URAs generate a single, detect an image. A typical example of commercial linear 
clear-cut peak, which corresponds directly to the location 45 CCD detector with several thousand pixels is the Kodak 
(displacement) of the image on the image detector. KLI-5001. Another example of a suitable commercial detec- 
tor is an EG&G Reticon RL4096, which comprises a linear 
OBJECTIVES AND SUMMARY OF THE array comprising a row of several thousand sequentially 
INVENTION addressable photodiodes. The sensor to which this invention 

..... r . . . , . , 50 is applied, and with which it is used, may be identical to 

* An objective of the present invention is to provide an ' , ... * .i. 

, ^ . . 1- 1.. f • . those that are known in the pnor art, or the sensor may be 

opucal apparatus which uses light from a point source to .... * . • *u ... • ti 

*^ . 1 -J 1 • -.1? 1- ui j^. • t. one that is not yet in the pnor art or not yet commercially 

create a nearly ideal unage with negligible distortion, such .... r^*^e / i . .i 

.1,- 1. , e ,u ' *u ' ' available, such as Imear CMOS (complementary metal 

that the lateral displacement of the image on the imaging . . ... \ u . j- j i. • \u d - ■ 

... ... -ii u J- *i I * oxide sihcon) photodiode arrays. It IS the means of bnnging 

pbotoscnsi .ve detector wUl be direcUy proportional to the radiation to the sensor and extraaing data from it, and 

tangent of the angle of the .ncoming hgbt rays. .^^^ ^^^^ constitutes this invention. 

Another objective of this invention is to provide an ^h^ ^ ^^^ „f ^ mwtniion comprises at least 

apparatus that creates an image that is m sharp focus ^^ of electromagnetic radiation and a receiver 

regardless of the distance of the pomt light source or lis ^f such radiation, where there is interposed between the 

angle with respect to the optical apparatus. ^ ^^^^ ^^^i^^, ^ p,^,^,i,y ^^^^^^ 

A further objective of this invention is to provide a linear preferably of various widths, where the slits are spaced apart 

image at least ten times the luminosity than that obtained in an irregular pattern. The slits should be wide enough and 

from a single simple slit while decreasing the relative efifects spaced far enough apart that diffractive effects are mini- 

of noise or artifacts of the system. mized. (Diffractive effects are undesirable, because the 

Another objective of this invention is to mitigate the 65 diffraction patterns are not stable with respect to changing 

effects of minor defects, dust, smudges, and fingerprints on source locations.) These slits may be literal openings 

the optics of the detector. through a self-supporting opaque material, or they may te 
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substantially transparent narrow elongated rectangular win- of ambient background illumination comprising a broad 

dows within a substantially opaque mask that is disposed on spectrum of wavelengths. In this embodiment, the filter can 

a transparent substrate. light from the point source is passed be the substrate for the aperture mask itself. Presiunably the 

through the sht pattern in such manner as to east a similar wavelength of the light source would itself be matched with 
pattern of undistorted Unear real images of the slits onto the 5 the wavelength at which the detector itself exhibits best 

detector. The axis of the detector is arranged to be substan- performance. For modern, typical CCDs this wavelength 

tially perpendicular to the lines of the image. The total might be in the range of the near-visible infrared wave- 

luminosity of the pattern is proportional to the cumulative lengths. Alternatively, there may be other, perhaps practical, 

width of the slits. While regularly spaced slits could be used, reasons for using a wavelength of light that is visible to 
they present a possible ambiguity when one or more of the lo human eyes. In any case, as new detectors are perfected, 

image lines at an end of the pattern are eclipsed. It is different wavelengths of hght may become more desirable, 

difficult, if not impossible to determine which end of the For larger scale applications, the techniques of this invention 

pattern is missing. Therefore, a non-repeating or pseudo- might employ other types of radiation, such as for example: 

random' pattern of slits is preferred. sound waves or radio waves, using appropriate slits with 

In a preferred embodiment of this invention there are at 15 aperture masks made of appropriate materials opaque to the 

least about ten (10) sUts, and the cumulative width of the lype of radiation used, 

slits is at least ten (10) times the width of a simple single slit ^^^^r.,^r„^., 

as has been used in the prior art. This apparatus exhibits BRIEF DESCRIPTION OF THE DRAWING 

distortion and depth of field comparable to those of a single xhe accompanying figures illustrate prior art as well as a 
slit, which are improved over that of a cylindrical lens 20 p^ferred embodiment of the present invention and thereby 

system. This improvement is particularly evident where the form a part of the specification. Together with the textual 

incident light is at more oblique angles. Yet, for certain description, these figures serve to explain the advantages 

designs, the total luminosity of the instant system could and principles of the invention. Each reference number in 

approach that of a cylindrical leas with baffies, the figures consistently refers to the same component of the 

Furthermore, there is an unexpected advantage of invented apparatus throughout all the figures as well as in the 

employing multiple narrow slits in comparison to either a textual description. 

cylindrical lens or a single sUt of any width. When precisely piG. 1 is a perspective schematic view of a lens imaging 

interpolating the fractional-pixel location of the image on sensor of the prior art. 

the detector, the most significant location information is 2 is a perspective schematic view of a slit imaging 

embodied m the edges of the image hne of each sht. sensor of the prior art. 

Therefore, a plurahty of slits contains multiples of the ^ . . - • ^ . . 

displacement information, much more redundant informa- ^ is a perspective schematic view of a sht imaging 

tion than is available from a simple image formed by a single according to this mvention. 

lens or a single slit. FIG. 4 is a plan schematic view of the operation of the 

In order to recover this extra image location information, P^^^^ °^ invention, 

much more computation is required than with either a lens FIG. 5 is a pattern of shts according to this invention, 

or a single slit. Nevertheless, this engineering trade-off will FIG. 6fl is a graph of the mathematical representation of 

be acceptable in general, because of the speed of electronic a pattern of slits according to this invention, 
computing processors, which is increasing each year, and ^ FIG. 6b is a graph of the mathematical representation of 

even at a decreasing cost. At present, the cost of a digital the corresponding reference pattern, 

signal processor generaUy is far less than the cost of a ^^^^ ^^ ^ ^^^^^^^ ^^^^^ ; 

precision multiple-lens optical system, and it is reasonable , . . , r . • r 

^^r^^r^f iUi^ Z.ct A;^^r.^i„ t« I^IG. 6d IS a graph of the correlation function as the 

to expect tills cost dispanty to continue. ...r,,^ , ^ 

^ t reference pattern is shifted left or nght from the peak 
Because the slit pattern forms a complex image on the 45 displacement. 

detector, it is preferred to use an imaging detector that has rTi- t • u- 1 * r j * . t. 

» 1 ^ u r J- . • t / u /^OT^ FIG. 7 IS a graphical presentation of data that show a 

a very large number of discrete pixels (such as a CCD or a . . , 1 £. j i_ . 

niiMric ««. A Tn, • 1 *i, * • *u typical correlation function and a continuous best-fit qua- 

CMOS photodjode array). The more pixels that are in the JJ^.- 

detector, the greater will be the accuracy of the resultant ninciion. 

determination, although at some point the price increases 50 ^ ^ S^'P^ corresponding to the graph of FIG. 6a, 

sharply. It is contemplated that deteaors with several thou- missmg apertures. 

sands of pixels will be suitable for this invention. A position 86 is a graph corresponding to the graph of FIG. 6d 

sensitive device (PSD), which only detects the location of showing the correlation of the image pattern of FIG, Sa with 

the cenlroid of the radiation impinging on it, is poorly suited reference pattern of FIG. 6b, 

to use in the practice of this invention, because it would not 55 FIG. 9 is a graph showing the height of the peak of the 

capture the redundant information coded into the complexity correlation function as the magnification (scale) of the image 

of the pattern. pattern changes. 

The electromagnetic radiation being discerned according FIG- 10 is a flowchart of the calculations associated with 

to this invention may be in the form of a broad spectrum of this invention. 

wavelengths, such as visible light. It need not be coherent or 60 FIGS. 1 through 4 include a number of common compo- 

monochromatic. Alternatively, narrow spectra wavelengths nents. One is a point source of light 10, illustrated by a tiny 

of radiation can be used to advantage. One means of infrared light-emitting diode (LED). Another common com- 

producing such a narrow radiation wavelength band is to ponent is a linear photosensitive image detector 14 that is 

employ an optical filter that passes only a narrow range of conventionally used to detect the real image (such as 31 in 
wavelengths which are matched to the wavelength of a 65 FIG. 1) of the point light source focused onto a row of many 

nearly monochromatic light source. Such embodiment of photosensitive pixels 16. The light source 10 generally will 

this invention generally can decrease the interference effects be at an oblique angle with respect to the image detector 14 
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and the optical oeaterliDe 11. In RGS. 1 through 3 the 
location of the light source 10 is depicted by a lateral offset 
distance 13 of the light source 10 from the optical oenterline 
11 and by an elevation distance 12 of the light source 10 
above the optical centerline 11. 5 

FIGS. 1 and 2 illustrate arrangements that represent the 
state of the prior art. FIG. 1 is a perspective view illustrating 
the conventional use of a cylindrical lens 21 or lens system 
to focus the light rays 15 from a point light source 10 into 
a linear real image 31, which is made to cross the photo- lo 
sensitive strip 16 of a linear image detector 14. FIG. 2 is a 
perspective view illustrating the use of a single straight slit 
22 to cast the light 15 from a point source 10 into a linear 
image 32 that is made to cross the strip 16 of a linear image 
detector 14. 15 

FIG. 3 is a perspective view of an arrangement of the 
apparatus of the present invention. It illustrates the use of an 
irregular pattern of multiple parallel slits 25 to convert the 
light 15 from a point source 10 into a pattern 35 of linear 
images that are made to cross the strip 16 of a linear image ^ 
detector 14. 

FIG. 4 is a plan view of an arrangement of the apparatus 
of the present invention that has been shown in perspective 
in FIG. 3. The light rays 15 are shown to flow through the 
irregularly spaced slits 25 to be converted to a corresponding 
congruent image pattern 35 on the photosensitive strip 16 of 
the linear image detector 14. 

FIG. 5 depicts an enlargement of a symmetrical, but 
otherwise irregular, pattern of slits, which is in fact an URA. 
In this figure, the black areas are intended to be the 
apertures, and the white areas opaque mask. (However, the 
negative of this could also work as a URA.) 

FIG. 6a is a graph of the mathematical representation of 
a specific pattern of apertures according to this invention. In 35 
this case, it is an URA based on the prime number N=97. A 
1 represents an aperture, and a 0 represents an opaque area. 
FIG. 6^? is a graph of the mathematical representation of the 
corresponding reference pattern. Here, a 1 represents where 
we expect to see image of radiation passing through the 40 
corresponding slit. A -1 represents where we expect to see 
the shadow of an opaque area. FIG. 6c is the graph of a 
reahstic, actual image pattern, which includes thermal noise 
and some diffraction artifacts, FIG. 6rf is a graph of the 
correlation between the image and reference patterns as the 45 
reference pattern is shifted left or right from the peak 
displacement. The horizontal scale of FIG. 6d is the dis- 
placement of the reference pattern (FIG. 6b) relative to the 
displacement of the actual image pattern on the detector. 

FIG. 7 is a graphical presentation of data that shows a 50 
typical mathematical correlation function 42 where the 
horizontal axis is the displacement (in pixels) of the image 
pattern 35 on the image detector 14, and the vertical axis is 
the correlation function value (the relative degree of match). 
It also shows a continuous best-fit quadratic function 44 with 55 
a peak 46, which estimates the sub-pixel displacement of the 
image pattern. The horizontal scale is the absolute displace- 
ment of the reference pattern in terms of the pixels of the 
image detector. 

FIG. 8fl is a graph corresponding to the graph of FIG. 60, 60 
but with missing apertures. This simulates a CCD with 
"dead" pixels or an URA with blocked apertures (because of 
dust or defects, for example), FIG. 86 is a graph correspond- 
ing to the graph of FIG. 6d showing the correlation of the 
defective aperture pattern of FIG. Sa with the reference 65 
pattern of FIG, 6b. This demonstrates the immunity of the 
URA to anomalous local defects. 



FIG. 9 is a graph showing the height of the peak of the 
correlation function as the size (scale) of the image pattern 
changes by a few percent, but the reference pattern remains 
the same size. 

FIG. 10 summarizes the calculations to compute angular 
location of the source with respect to each sensor and to 
compute the XYZ location of the source assuming that there 
arc a sufficient number of sensors at appropriate locations. 

The parts in the various figures are numbered throughout 
as follows: 

10 point source of light 

11 optical centerhne 

12 elevation distance of light source above centerline 

13 lateral oflfeet distance of light source from centerhne 

14 photosensitive image detector (photodetector) 

15 ray(s) of light as they pass from the point source to the 
photodetector 

16 row of photosensitive pixels 

27 cylindrical lens 
22 straight slit 

25 irregular sht pattern 

28 multiple parallel irregularly -spaced slits 

31 real image of a point source through cylindrical lens 

32 real image of a point source through a single slit 

35 real image of a point source through irregular slit 
pattern 

42 correlation between acmal and reference images for 

various displacements 
44 best-fit interpolation function (such as a quadratic 

polynomial) 
46 peak value of best-fit interpolation function 

DETAILED DESCRIPTION OF A SPECIFIC 
EMBODIMENT 

The invention will be described below with reference to 
the Figures and the numbered individual components 
therein. In the description below, the specific construction, 
the count, and the arrangement, of the components are 
intended for clarity of illustration. Other, and perhaps even 
belter, arrangements of the components constitute alterna- 
tive specific embodiments of the same method and appara- 
tus. FIGS. 3 and 4 illustrate this invention in the use of an 
irregular pattern 25 of multiple straight and parallel slits 28 
to convert the light from a substantial point source 10 into 
an image pattern 35, which is caused to cross the photosen- 
sitive image detector 14. This invention operates somewhat 
differently than prior art, in that the prior art used a con- 
verging cylindrical lens. Instead of gathering and focusing 
more light intensity onto the same fine image location on the 
photosensitive strip 16, the instant invention generates mul- 
tiple image lines from multiple slits 28. Thereby the total 
amount of light impinging upon the image detector is many 
times the amount of hght available through the use of a 
single slit. To reduce diffraction effects, these slits 28 are 
large enough (suitably about 0.3 to 2.0 mm in various 
widths) and separated by large enough opaque areas (also 
about 03 to 2.0 mm). The photosensitive strip in this case is 
a row of several thousand CCD pixels. Alternatively, these 
pixels can be those of a comparable linear CMOS image 
detector. The large number of pixels allows the complexity 
and details of the image pattern to be detected. The image is 
then processed sequentially, pixel by pixel, using a conven- 
tional digital signal processor (DSP) to read, identify, and 
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locate the displacement (of the center) of the image pattern 
in the midst of electronic noise using standard signal pro- 
cessing techniqties. 

All the slits should have substantially parallel sides. The 
straighter and more parallel the sides are, the better will be 5 
the definition of the pattern 35 and the more accurate will be 
the resultant determination of the centroid of the point light 
source. These slits, however, are of irregular width and 
spacing. The reason for the irregularity is to restrict or 
prevent ambiguity of the displacement of the pattern when jq 
the images from some of thie individual slits are missing 
from one or both ends of the image pattern. In fact, ideally 
designed irregularity will allow part of the image to extend 
beyond the end of the photosensitive strip but will still allow 
the displacement of the image's center to be unambiguously 
determined. The DSP, using the actual image pattern, the 
reference image pattern, and conventional correlation 
calculations, in effect sums the intensities of the multiple 
images while random errors, due to background noise for 
example, tend to cancel statistically. In particular, the total 
image signal amplitude grows faster than the average total 
noise so that the effective signal-to-noise ratio (SRN) 
increases with the square root of the total amplitude. 

The brute -force way to determine the location of the 
image, that is, the displacement of its center relative to the 25 
detector, is to correlate the actual image pattern with a 
reference pattern. The correlation can be performed by 
incrementally displacing the reference pattern from left to 
right by each possible pixel displacement until the correla- 
tion maximum (best match) is found. The displacement with 30 
the best match (the peak of the correlation function) then 
corresponds to the tangent of the angle which the point 
source forms with respect to the sensor's normal vector. This 
is illustrated in FIG. 6d, which plots the correlation between 
an ideal image pattern and the reference pattern for various 35 
displacements of the reference pattern relative to the image 
pattern. The correlarion computation is a very well- 
established and common technique in digital signal process- 
ing and image processing and is described in almost any 
textbook on the subjects. 40 

The correlation function provides values only at discrete, 
whole -pixel displacements. See FIGS, 6d and 7 for 
examples of the plot of such discrete correlation calcula- 
tions. In FIG. 7 the stair-step graph 42 of the correlation 
function shows the effect of displacing the reference pattern 45 
by different amounts as measured in whole pixels. For 
example, with the reference pattern shifted by N pixels, the 
correlation function yields the highest value. A shift of N+l 
pixels yields a lesser value. The precise actual displacement 
of the image pattern, as measured to a fractional pixel 50 
(sub-pixel), can be estimated by using well-known best-fit 
polynomial interpolation techniques. This can be done by 
using the pixel displacement with the maximum correlation 
value and by using several of its immediately neighboring 
correlation values on each side. The dotted line 44 shows 55 
how the peak 46 of a best-fit continuous quadratic polyno- 
mial might be used to interpolate the exact displacement of 
the image in terms of fractions of a pixel (for example, 
N+0.41 in FIG. 7). Other techniques may be used to inter- 
polate the actual displacement, such as computing the cen- 60 
troid of the graph 42. 

The interpolated displacement of the center of the image, 
in effect, uniquely determines a plane which contains the 
point source of li^t 10, the centerline of the aperture pattern 
25 and parallel to the apertures 28, and the center of the 65 
image pattern 35, The location of the point radiation source 
is the intersection point in 3-D space of all the planes from 
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all the multiple sensors. If there are more than three planes 
(that is, more than three sensors), the invention would use a 
"best-fit" method to determine that point (such as the 
common method of finding the point which minimizes the 
sum of the squares of its distances to all the planes). This is 
because the measurement will have some small error, and 
fiirther because the point somx^ of the radiation is really not 
a point source in the strictest sense of the term. Rather it is 
a close approximation of a point radiation source and 
therefore the muhiple planes computed from the multiple 
sensor displacements may not intersect in exactly one point 
in space. Given enough sensors at appropriate known loca- 
tions and orientations, the coordinates of this point in some 
fixed coordinate system are unique and may be computed. 
The mathematical details for one way to perform this 
calculation arc presented in the Burton and Sutherland paper 
cited previously. 

In practice a digital signal processor (DSP) associated 
with each image detector would compute the correlation 
function to determine the exact image displacement. Then a 
standard floating point microprocessor, such as an Intel 
Pentium processor, would compute the XYZ coordinates of 
the location of the point source in the 3-D spatial coordinate 
system. The latter processor would be provided with the 
interpolated digital image displacement value from each 
sensor's DSP as well as calibration information such as the 
locations and orientations of the sensors 

While almost any random irregular pattern of slits would 
be adequately effective in the practice of this invention, 
some patterns have better correlation characteristics than 
others. Ideally, the correlation computation will return a zero 
value for all displacements of the reference image except 
those near the correct displacement. See FIG. 6d as an 
example. Near the correct displacement (say at pixel N), the 
correlation value should increase sharply to a peak 46, but 
not too sharply, since it is desirable to use some interpolation 
fixnction in order to estimate the sub-pixel location of the 
peak 46 as shown in FIG. 7. The least -squares best-fit 
quadratic would be such an interpolation function. 

The pattern need not consist simply of transparent aper- 
tures in an opaque mask as shown in FIG. 3 through FIG. 5. 
Alternatively, it could consist of a pattern of strips of varying 
transparencies or even a wave -like pattern of a continuously 
varying transparency. Such continuously varying patterns 
will render better interpolation if the pixel size approaches 
or exceeds the size of the aperture or strip width. This is not 
the case for the preferred embodiment described herein, 
where each slit or intervening opaque area casts an image 
many pixels wide. High contrast patterns arc easier to 
generate and less costly to manufacture in quantity with 
good precision and quantity. 

As noted above, the field of x-ray imaging uses a related 
concept of a pattern of multiple pinhole apertures to increase 
total film exposure while maintaining sharp images. Such 
coded apertures traditionally use irregular 2-dimensiooal 
patterns of round or rectangular holes. The resulting over- 
lapping images are used for the purpose of creating complete 
images of complex x-ray illuminated objects rather than for 
the purpose of precise, sub-pixel displacement measurement 
to localize a point source of radiation. 

A uniformly redundant aperture (URA) pattern is a spe- 
cific kind of coded aperture which possesses the desired 
characteristics described previously. Furthermore, it has an 
average density of about 50%, so it will pass much light to 
the detector. Therefore, the preferred embodiment of the 
present invention employs an URA-type pattern, but in only 
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one (1) dimension. An example of a uniformly redundant centroid. Then, instead of the FFT computation above, the 

aperture pattern for our 1-dimensional purposes is shown in correlation function could be computed directly at only 

FIG. 5. about a dozen displacements and the interpolation computed 

A URA pattern (for the instant 1-dimensional purpose) is correlation values. In this case, the image 

based on some minimal unit slit width W (such as 03 mm, 5 displacement could be evaluated m a tunc interval that is 

r ,N ... Li -nw; -J L equal to or less than that needed for the FFT approach, 

for example), and the whole pattern is P*W mm wide where n •* • j *u . j- . e .a. 

n ' ' , L j.L J L Normally it is assumed that the distance from the pomt 

P IS a prime number. Each sht width and each opaque ^^^^^^^ ^^^^^ ^ separ^Uion 

separation between two sbts is some multiple of W. If P ^^^^^^ ^^^^^ '^^^^^^^^ ^ 

umt-widc contiguous rectangular ^^^^^^^^ ^^^^ ^y^^^ the aperture mask, the 

-(P-l)/2 on one end, through 0 m the middle, and up to lo jmage pattern on the detector grows somewhat in size. This 

(P-l)/2 at the other end, then strip S is transparent if and must be accounted for, because a growth of just a few 

o°^y ^ percent in the image significantly degrades the height and 

width of the peak in the correlation function. This is depicted 

I -s modulo P, for some integer L ^ t^c graph of FIG. 9 for the case of a 15 mm wide URA 

T^,<.,^f.^.^ n>^i\n ct.;«c .k* with P-97. One solution to this problem is to scale down 

If p f ^ t ansp^ent, and the est are ^^^^^j ^ ^^^^ ^ 

opaque. If P=4*I+1 for some mteger I, then the pat ern is ^^ich is ^he same as the width of the aperture mask itself, 

symmetrical. Otherwise the pattern is anti-symmetrica^. That ^ elementary way of doing this is simply to delete M 

is, each transparent stnp on one side of the centerhne is ^.q^^Uy separated pixels from the image, where M is the 

mirrored by an opaque stnp on the other side of the 20 difference in sizes of the actual image and of the aperture 

centerline. Groups of contiguous transparent strips form mask (as measured in pixels), until the total remaining image 

composite apertures of various widths. Groups of contigu- is of the correct size, A more precise solution is to perform 

ous opaque strips produce varying separations between the correlation calculations with reference patterns of several 

apertures. The best choice of P and the choice of the exact sizes and then use the resulting correlation function with the 

width of the strips depends on the length of the detector, the 25 highest peak (that is, the best match). Note that this would 

number of pixels, the desired field of view, the distance yield an estimation of the distance of the point source from 

between the aperture mask and the detector, and diffraction the sensor. That estimate could be used to verify the final 

limitations for the chosen wavelength of light. For example, calculation of the location of the point source, 

if the detector has 4096 pixels each of which is 8 microns Another possible solution to the problem of a magnified 

square (totaUng 33,77 mm long) and it is desired that the 30 ™age is to recompute (the FFT of) the description of the 

image pattern cover about 1/4 of the detector (or 1024 pixels), reference pattern dynamicaUy, based on the width of the 

then the aperture pattern should be 8 mm wide. For Po23, p^i^m. This however assumes that the whole 

the strips should be about 0.348 mm wide. For P=17, the is present so that its width can be determined directly, 

strips should be about 0.471 mm wide. For purposes of this ?i'w°'!^f ''''^Y^^^''^^}^'' reference pattern vndih 

* • 1 ' r J L • P*W could be determmed by estimating the Width W of the 

mvention, a symmctncal pattern is preferred because it is ,e * • ^ .1. v. ■ i.- ^ • -1.1 ^1 t 

, ' ... J . umt stnp from the slit images which are visible. Clearly, 

more hkely to produce consistent accuracy under varying other variations of these optiSiized methods are possible and 

conditioi^. Tins is because there IS no possibility of intro- invention, 

ducing a bias toward one side or the other m the presence of Two-dimensional aperture arrays and area image detec- 

diffraction artifacts. the methods described above. For 

When computing the correlation function as shown in 40 example, the referenced paper by E. E. Fenimore applies 

FIG. 6d, thousands of multiplications and additions must be 2-D coded apertures to 2-D imaging such as with x-rays, 

done for each of thousands of discrete displacements of the However, the computational cost of doing the general 2-D 

reference pattern. This results in a total of milUons of correlation computation (with or without FFTs) prohibits 

arithmetic operations. Therefore, some optimized algorithm 2-D displacement calculation in real-time with current 

or short cut is desirable. One preferred optimization is to 45 technology, since there are at least 250 thousand pixels on a 

perform a fast Fourier transform (FFT) on the actual image typical two-dimensional area CCD. Thus there are substan- 

and an FFT on the reference pattern. The latter FFT need tially more pixels to be processed and there are substantially 

only be computed one time, since the reference pattern is displacements of the reference image at which to 

constant. After multiplying the two FFTs element by compute the correlation. Of course, given an estimate of the 

element, the inverse FFT of the product is computed, the 50 2-D displacement of the image on the detector, one can use 

pixel displacement with the peak value is found, and finally short cuts as in the 1-D case to drastically reduce the 

the results are interpolated using several neaiby values. This computation time. In the 2-D detector case, the need for this 

computation will take only tens of thousands of arithmetic ^inj of alternative optics is not as important, because 

operations: standard 2-D spherical lens systems are more practical in 

^PP^i.pj^,p.,p.p.^^ . 55 that case. Standard lens systems are less problematic, are 

readily available in large varieties, and are more highly 

where developed than cylindrical lenses. 

C is the correlation value at each integral pixel Although this invention has been described above with 

displacement, * reference to a preferred embodiment, anyone skilled in the 

. . . , ... ^ ^ 60 art can readily visualize alternative embodiments of this 

D IS image pixel amphtudes on the deteaor, ^^^^^^^^ Therefore, the scope and content of this invention 

E is +1 for an aperture in the reference pattern and -1 are not limited by the foregoing description. Rather, the 

elsewhere, and scope and content are to be defined by the following clauns: 

* is element-by-element complex multiplication What is claimed is: 

In the case where the actual image pattern is known to be 65 1. In an apparatus for determining a location of at least 

whole, its approximate displacement could be estimated one substantially point-like radiator of energy in a three 

from the pixel locations of its extremities or from its dimensional volume comprising; 
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said at least one substantially point-like radiator of 
energy; 

at least one detector adapted to be disposed in a known 
location within said three dimensional volume, and 
adapted to be positioned to intercept energy radiating ^ 
from said substantially point-like radiator, wherein said 
detector comprises at least one array of elements 
adapted to detect said radiation and to convert said 
delected radiation into output signal^ and 

at least one pattern member, adapted to be disposed 
between said array and said radiator, comprising plural 
gradations of transparency, wherein said pattern is 
adapted to spatially modulate radiation that is incident 
thereon as a frinction of said gradations in transparency, 
and lo project said modulated energy onto said array in 
the form of at least one real image; and 

a signal processor adapted to convert output signals 
representative of said at least one modulated real image 
into a direction of said radiator relative to an assembly 
of said detector and said pattern member; 
the improvement that comprises: 

at least one substantially ideal reference image of said 
pattern member on said array; 

said signal processor being adapted to compare said 25 
modulated image with said reference image; 

said signal processor comprising a discrete mathematical 
correlation function that is adapted lo characterize the 
degree of match between said modulated image and 
said reference image for a plurality of discrete displace- 30 
menls of said reference image with respect to said 
array, 

said signal processor being adapted to identify at least one 
of said displacements that produces a correlation that is 
characterized by a well defined maximum value; and 35 
said signal processor being adapted to convert said 
displacement(s) producing a maximum value into a 
direction of said radiator in said three dimensional 
volume with respect to the known location of said 
detector in said three dimensional volume. 40 

2. The improved apparatus of claim 1 further comprising 
said signal processor being adapted to estimate the displace- 
ment of the maximum value of a continuous function, that 
equals the correlation function at said displacements inter- 
polated between at least two discrete displacements of said 45 
correlation function; and 

said signal processor being adapted to convert said esti- 
mate into an interpolated, and therefore more precise, 
direction of said radiator in said three dimensional 
volume. 50 

3. The improved apparatus of claim 1 wherein said plural 
gradations of transparency comprise a plurality of slits of 
varying width in an otherwise substantially opaque field. 

4. The improved apparatus of claim 3 wherein said array 

is substantially one dimensional and comprises a linear strip 55 
of said radiation converters, wherein long dimensions of 
said slits are substantially normal to said linear strip. 

5. The improved apparatus of claim 4 comprising a 
plurality of said arrays at least two of which are disposed 
substantially normal to each other; and said signal processor 60 
is adapted to determine directions of said radiator with 
respect to said arrays in at least two dimensions correspond- 
ing to the dispositions of said normal arrays. 

6. The improved apparatus of claim 4 comprising a 
plurality of said arrays at least three of which are disposed 65 
not all parallel to each other and said signal processor is 
adapted to determine directions of said radiator with respect 



to said arrays in at least three dimensions corresponding to 
the dispositions of said arrays whereby enabling said signal 
processor to determine the substantial point location of said 
radiator in said three dimensional volume. 

7. Tlic improved apparatxis of claim 3 wherein the width 
of said slits are not digitally related to each other. 

8. The improved apparatus of claim 3 wherein the widths 
of said sUts correspond to a uniformly redundant array 
(URA) pattern. 

9. The improved apparatus of claim 3 wherein said slits 
are of random width or random spacing. 

10. The improved apparatus of claim 3 wherein said slits 
are disposed symmetrically about an arbitrary line that is 
contained in said real image and is transverse to said array. 

11. The improved apparatus of claim 1 wherein said 
radiation is electromagnetic radiation. 

12. The improved apparatus of claim 11 wherein said 
radiation is infra red. 

13. The improved apparatus of claim 1 wherein said array 
comprises a plurality of photodc lector elements. 

14. The improved apparatus of claim 13 wherein at least 
some of said elements are linearly aligned. 

15. The improved apparatus of claim 14 wherein said 
array comprises a charge coupled device (CCD). 

16. The improved apparatus of claim 15 wherein said 
CCD comprises at least one linear row of photosensitive 
pixels. 

17. In a method of determining the direction of at least one 
unknown point in a three dimensional volume, relative to a 
different, known location in said volume, which comprises: 

disposing at least one substantially point source energy 
radiator in said three dimensional volume in known 
spatial relationship to at least one of said point(s); 

disposing at least one detector, that is adapted to detect 
energy radiated by said radiator and impinged upon at 
least one linear array of elements that are adapted to 
detect said radiation and lo convert detected radiation 
into output signals that are proportional to strength of 
said impinged energy on the elements of said array, 
wherein said detector is in known spatial relationship to 
said three dimensional volume; 

operatively associating said detector(s) with a signal 
processor that is adapted to convert said output signals 
to a direction of incidence of said radiation on said 
detector(s); 

the improvement that comprises: 

interposing a pattern member, comprising plural areas 
having gradations of transparency to said radiation, 
between said radiator and said array; 

radiating energy from said radiator to said array through 
said pattern member such that a modulated image of 
said pattern member is disposed on said array; 

generating an output of said array that is representative of 
at least one displacement of a substantially ideal ref- 
erence image of said pattern member; 

in said signal processor, comparing said modulated image 
with at least one of said displaced reference images 
until a best match is found between a specific displaced 
reference image and said modulated image; and 

in said signal processor, from said best match, determin- 
ing a direction of said radiator with respect to said 
detector. 

18. An improved method as claimed in claim 17 further 
comprising repeating said direction determination with said 
detectors at different spatial relationships to said volume, 
sufi&cient to establish the three dimensional coordinates of 
said unknown point. 
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19. An improved method as claimed in claim 18 wherein 
said array is a plurality of linear aligned plural elements that 
are adapted to convert impinged energy into output signals 
and fiirther comprising repeating said improved method with 

at least three arrays disposed not all parallel to each other. 5 

20. An improved method as claimed on claim 18 wherein 
said plural gradations of transparency comprise a plurahty of 



18 

slits of differing, non-digital widths that are longitudinally 
disposed normal to said array; and fiirther comprising caus- 
ing said modulated image to comprise a set of lines trans- 
verse to said array. 

4* « ♦ 4* * 
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